Three dimensional resonating valence bond liquids and their excitations 
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We show that there are two types of RVB liquid phases present in three-dimensional quantum 
dimer models, corresponding to the deconfining phases of 17(1) and Z2 gauge theories in d = 3 + 1. 
The former is found on the bipartite cubic lattice and is the generalization of the critical point 
in the square lattice quantum dimer model found originally by Rokhsar and Kivelson. The latter 
exists on the non-bipartite face-centred cubic lattice and generalizes the RVB phase found earlier 
by us on the triangular lattice. We discuss the excitation spectrum and the nature of the ordering 
in both cases. Both phases exhibit gapped spinons. In the U(l) case we find a collective, linearly 
dispersing, transverse excitation, which is the photon of the low energy Maxwell Lagrangian and we 
identify the ordering as quantum order in Wen's sense. In the Z2 case all collective excitations are 
gapped and, as in d = 2, the low energy description of this topologically ordered state is the purely 
topological BF action. As a byproduct of this analysis, we unearth a further gapless excitation, the 
piOn, in the square lattice quantum dimer model at its critical point. 



I. INTRODUCTION 

Fractionaliscd phases, i.e. phases where low energy 
excitations exhibit fractional quantum numbers, have 
played an important role in condensed matter physics 
over the past two decades. Probably the two most salient 
examples are provided by the quantum Hall effect, where 
the charge of a quasiparticles can be a rational fraction 
of the fundamental electronic charge^ and by Anderson's 
suggestion of spin-charge separation being at the root of 
high-temperature superconductivity^ 

This latter proposal has spawned a set of theories 
collectively known as resonating valence bond (RVB) 
theories. 3,4 In its short-range incarnation^ RVB theory 
uses dimers as its basic degree of freedom, which are 
cartoons for a pair of adjacent electrons forming a sin- 
glet ('valence') bond. At zero doping, when all sites are 
singly occupied, all electrons are involved in a singlet 
bond with a neighbour. The Hilbert space of this theory 
is thus given by the classical dimer coverings of the lat- 
tice. Resonance moves between dimer coverings provide 
a quantum dynamics for the dimers. 

The original hope was that the resulting quantum 
dimer model would have a liquid phase, the RVB liquid, 
hypothesized as an alternative to Nccl order in antiferro- 
magnets by Anderson and Fazekas in the early seventies. 5 
This state would exhibit fractionalized 5=1/2 spinons 
and as a consequence holes doped into it would undergo 
spin-charge separation and the resulting charged spinless 
quasiparticles ('holons') would Bose condense to form a 
superconductor. 

However, it turns out that the square lattice quan- 
tum dimer model appropriate for the cuprates ex- 
hibits only solid phases with the exception of one 
critical point 4&L£i2i These properties are now well- 
understood within the framework provided by the height 
representation!!! for dimer coverings^!*!^ By contrast, for 
dimer models on non-bipartite lattices, the present au- 
thors have demonstrated that an RVB liquid phase is 



possible*! 3 - The different outcomes on bipartite and non- 
bipartite lattices can be rationalised, following the ideas 
of Fradkin and Shenker on lattice gauge theories^* via 
the absence (presence) of a deconfmed phase in {7(1 ) (Z2) 
gauge theories in 2+1 dimensions i^ik We note that the 
RVB phase on non-bipartite lattices is best character- 
ized as a topologically ordered phase!!! whose low energy 
theory is the purely topological BF theory>!£ 

This understanding of the situation in 2+1 dimensions 
has led us, in the present work, to investigate the be- 
haviour of quantum dimer models in d = 3+1, specifically 
on the bipartite cubic lattice and the non-bipartite face- 
centred cubic lattice. We find that both models exhibit 
RVB phases with deconfmed spinons, which correspond 
to the Coulomb phase of the U(l) gauge theory (ordi- 
nary electromagnctism) and the deconfined phase of the 
Z 2 gauge theory, respectively. The FCC lattice RVB liq- 
uid is topologically ordered. It has a gap to all excitations 
and its ground state degeneracies and topological interac- 
tions between spinons and vortex loops are encoded in the 
topological 3+1 dimensional BF action. By contrast, the 
cubic lattice RVB phase supports a gapless, transverse 
collective mode. This can be identified as the photon of 
the low energy theory, which is now a non-topological 
Maxwell theory. Consequently, we identify the ordering 
of this phase as quantum order in Wen's sensed reserv- 
ing the term topological order for cases where the low 
energy theory is purely topological. As a byproduct of 
this investigation, we have also revisited the theory of the 
critical point of the square lattice dimer model where we 
find, in addition to the resonons which are the analogs 
of our photons in d = 2, a further gapless excitation, the 
piOn, which signals the incipient crystalline order. We 
note that our main results were announced in Ref . Il9i 

From the foregoing it follows that a deconfined phase 
in a valence bond dominated regime is easier to achieve 
in d = 3 than in d = 2, where its existence is ruled out 
for bipartite lattices. This encouraging fact is, however, 
counterbalanced by the difficulty of stabilising a valence- 
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bond dominated phase in comparison to a Neel phase, 
which becomes increasingly hard as the coordination of 
the lattice grows. There is also the possibility of real- 
izing dimer models in other contexts - e.g. the Santa 
Barbara group has shown that multiple dimer models 
arise when frustrated Ising magnets are imbued with a 
ring-exchange dynamics. Indeed in Ref.|23the reader can 
find a parallel discussion of a Colulomb phase in multi- 
dimer models on the diamond and cubic lattices. Alter- 
natively, dimer models can arise as Bose Mott insulators, 
electronic Mott insulators at fractional fillingsSi and in 
mixed- valence systems on frustrated lattices^ 

We would be remiss if we did not note that quite gener- 
ically U(l) gauge theories of Heisenberg magnets should 
be expected to predict a Coulomb phase in d = 3. Specifi- 
cally, a treatment of three-dimensional quantum magnets 
in the bosonic large- N theorySiS4 at a fixed (and small) 
number of bosons per flavour, supplemented by l/N cor- 
rections yields results along the lines presented here for 
dimer models. The dimer models themselves can be ob- 
tained, at order l/N, in the large- N limit taken with 
a fixed boson number, and hence a vanishing number of 
bosons per flavor^ It is a remarkable fact that these two 
very different limits yield the same physics. The equally 
remarkable fact that they also agree in d = 2 was estab- 
lished previously. It appears safe to conjecture that this 
agreement will hold in higher dimensions as well. 

Finally, we note that Wen has recently urged that the 
physical photon be viewed as the consequence of a quan- 
tum ordered vacuum^ The models he has considered 
that give rise to an emergent photon have considerable re- 
semblence to the cubic dimer model considered in this pa- 
per. Indeed, dimers and monomers have a natural string 
interpretation when superposed on a reference configura- 
tion. 

In the following, we first briefly summarise the known 
properties of the quantum dimer model which are of rel- 
evance to us. We then describe the RVB liquids, first on 
the simple and second on the face-centred cubic lattice. 
In both cases, we discuss in turn correlations, topological 
properties and the excitation spectrum. 



II. QUANTUM DIMER MODELS: 
GENERALITIES 

On a general lattice, allowed moves between two differ- 
ent hardcore dimer configurations consist of interchang- 
ing alternately occupied and empty links forming a closed 
loop, known as a resonance loop. If we denote the 
presence (absence) of a dimer by an Ising link variable, 
(7 X = +1( — 1), such a resonance move is generated by the 
kinetic operator 

fni/2 \ 

f ° = n °M-i°i + k - c - > (2-i) 
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FIG. 1: Master phase diagram for Rokhsar Kivelon quantum 
dimer models. Special points are highlighted: the Rokhsar- 
Kivelson (RK) point, the transverse field (TF) point, and the 
maximally flippable (MF) point. 

where the loop, denoted pictorially by o, contains n; 
links, and are the raising and lowering operators for 
to <j x . A corresponding projection operator, 

Po = T D 2 (2.2) 

determines whether dimers are arranged on a given loop 
so as to be able to carry out a resonance move; such loops 
we call flippable loops. 

The quantum dimer Hamiltonian proposed by Rokhsar 
and Kivelson considers only the shortest possible reso- 
nance loops. It contains two terms, each involving one 
of these operators. The first carries out resonance moves 
with a kinetic energy t , and the second exacts an energy 
cost w for a flippable loop»2£ 

Hqdm = ^2-tof + VoVo (2.3) 

o 

Some features are common to the phase diagrams of 
quantum dimer models on most lattices. We have de- 
picted a generic phase diagram in Fig. ^ For v/t > 1, 
it is favourable not to have any possible flippable loops; 
configurations satisfying this requirements can usually be 
found, and the resulting phase is called the staggered 
phase, a name borrowed from the appearance of one of 
the inert square lattice configurations. For v/t = — oo, 
by contrast, one wants to maximise the number of flip- 
pable loops. In both cases, there are no quantum flu- 
cutations. These phases are confining, and confinement 
typically exists for most values of v/t, with the possibility 
of transitions between different confining phases. 

From the point of view of deconfined, liquid^ phases, 
the region v/t < 1 is the most interesting, because the 
balance between kinetic term and potential term overly 
favours neither high or low densities of resonance loops. 
The ground state wavefunction can then be spread mod- 
erately uniformly over most of the classical configuration 
space; if the ensemble of classical configurations is suffi- 
ciently disordered, this can then lead to an RVB liquid 
phase. 

This statement can be made precise at the Rokhsar- 
Kivelson (RK) point v/t = 1. Define a sector as the set 
of dimer coverings connected by T . Then the equal am- 
plitude superpositions of dimer coverings in each sector 
is a ground state at the RK point. Hence equal time 
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correlations can be computed as correlations of the uni- 
form fugacity classical dimer model on the same latticed 
These can often be obtained analytically - exactly (es- 
pecially in d < 2) or at least asymptotically - or from 
Monte Carlo simulations. 

Consequently, advances in determining classical corre- 
lations give valuable information on the quantum model. 
Indeed, it was the solution of the classical dimer model 
on the triangular lattice with its short-ranged corre- 
lations that allowed the triangular RVB liquid to be 
discovered^ 

Recently, Huse, Krauth and the present authors have 
discussed the classical dimer correlations on the cubic 
and FCC lattices^ In the next section we will build on 
this work to establish the presence and nature of RVB 
liquid phases on these two lattices. 



III. CUBIC LATTICE: AN RVB U(l) LIQUID 

The general structure of the phase diagram on the cu- 
bic lattice is indeed as sketched in FigQ] For v/t > 1, the 
ground states are the non-flippable ('staggered') configu- 
rations, such as the ones obtained by appropriately stack- 
ing the two-dimensional staggered configurations. For 
v/t large and negative, there are six maximally flippable 
columnar configurations, again obtained from stacking 
planar columnar configurations. Apart from these two 
crystalline phases we have been able to establish the ex- 
istence of a liquid, RVB phase for v/t < 1. This is a 
"Coulomb" phase, characterised by algebraically decay- 
ing dimer ( "magnetic field" ) correlations and whose spec- 
trum contains a linearly dispersing transverse collective 
excitation ( "photon" ) , a gapped topological defect ( "elec- 
tric monopole" ) in addition to Coulombically interacting 
spinons ("magnetic charges"). We now sketch the deriva- 
tion of these results. 



A. RK Point 

We begin with the properties of the RK point, v/t = 1. 
The classical dimer model on the cubic lattice was studied 
m Ref.H It was shown there that a useful parametriza- 
tion of the dimer configurations is in the language of 
solenoidal magnetic fields. Briefly, on the bipartite cubic 
lattice, on can think of a dimer as a magnetic flux, B 
of strength 5/6 = 1 — 1/z pointing from sublattice A to 
sublattice B; here, z is the coordination of the lattice. 
Identifying an unoccupied bond with flux —1/6 = —1/z, 
one has V ■ B = - magnetic charges are excluded. 

Upon local coarse-graining, analogous to what is 
done for models with height representations in two- 
dimensions, long wavelength configurations acquire an 
entropic weight quadratic in B and this allows compu- 
tation of the long distance correlations of the classical 



model. These are purely dipolar, 

(B^BjiO)) = l 3XiXj ~ r25ij (3.1) 

where K can be determined via Monte Carlo 
simulations M The language of magnetic fields also al- 
lows a characterization of the sectors. For a cube with 
periodic boundary conditions, the flux through any sur- 
face that wraps around it is invariant under local dimer 
moves and is invariant under lattice translations of the 
surface. In particular, if we let £i be planes perpendic- 
ular to the cubic unit vectors e^, the fluxes $f through 
them are the invariants that characterize a given sector of 
the dimer model. This flux is also known as a U(l) wind- 
ing number^ For an L 3 cube these range in magnitude 
from to L 2 /2. The form of the correlations given above 
is the average over all sectors, but also holds individually 
for all sectors with $ s = |E,*f e i| ~ O(VZ), which 
dominate the sum at large L. Finally, we note that the 
"staggered" configurations maximize $ B , or more accu- 
rately, the average magnetic field strength. 

At the RK point equal amplitude superpositions in 
all magnetic flux sectors constitute degenerate ground 
states. As we will see below, for v/t < 1 this degeneracy 
is lifted and only the sectors with $ B ~ o{s/L) will be 
important in the low energy analysis — so we will not need 
more classical information than given in Eq. 13.11 above. 
For this reason it is sufficient to focus on the zero flux 
sector at the RK point for now. 

As the exact ground state wavefunction is known, we 
can use the single mode approximation (SMA) for dimer 
models, pioneered by RK in their original paper— to get 
a good description of the excitation spectrum. We re- 
mind the reader that the SMA will provide a rigorous 
upper bound on the excited state energy at each momen- 
tum. As RK's explanation was rather compressed and 
their treatment incomplete in one important respect in 
the d = 2 case, we will provide a fuller account between 
this section and the appendices. 

Single Mode Approximation: Let |0) denote the 
ground state of the dimer model under consideration. As 
above, we define cr|(r) as the Pauli spin operator the 
eigenvalues ±1 of which correspond to presence and ab- 
sence of a dimer on the link at location r. Here, we have 
added the 'direction' vector f for clarity to distinguish 
between dimers pointing in the different possible direc- 
tions. 

Next, we define the Fourier transform of the dimer 
density operator, 

a?(q) = £af(r)e*p(iq.r). (3.2) 

r 

Acting with this operator on the ground state yields a 
state 

|q,f)=Sf(q)|0) (3.3) 
with nonzero momentum q so that 

(0|q,f)=0. (3.4) 
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Provided that 



(q,f |q,f) ^0 



(3.5) 



this excitation has a variational energy of 



p., n ^ <- (0| [^(-g), [H QDM ^m]] 10) _ /(<& ^ 
<0| of (— q)of (q) |0> a(q) 

where the numerator and denominator in the final line 
are the oscillator strength - /(q) - and structure factor - 
s(q). These can be evaluated as ground state expectation 
values. 

If it now so happens that the density <7?(qo) is a con- 
served quantity, then [Hqdm , 5?(qo)] = 0, and thence 
£J(qo,f) = 0. From the orthogonality relation (Eq. 13. 4|) 
it then follows that this is a variational upper bound on 
the energy of the 'excited' state |qo,r). Similarly, the 
behaviour of /(q) near qo can be used to determine a 
bound on the dispersion of the soft excitations. Further, 
and again only if Eg. 13.51 holds, a smooth behaviour of 
/(q) while s(q) diverges, can also be used to infer gap- 
less excitations. Indeed, this is the classic signature of 
incipient order. 

RK identified a conserved quantity for the square lat- 
tice, which generalises to the cubic case as the density 
of dimers pointing in a given direction (f = x, say) at 
wavevector qo = (q x ,Tr,ir). This follows from the fact 
that the quantum dynamics always creates and destroys 
pairs of neighbouring dimers. This implies that the oscil- 
lator strength, /(qo), vanishes. In fact, at a wavevector 
qo + k, one finds 



/(k) cx (k x f) 2 



(3.7) 



as outlined in Appendix^ This form holds for all values 
of — oo < v/t < 1. 

This, however, does not lead to lines of zero energy 
excitations because the structure factor also vanishes for 
all qo with q x ^ ir - Ea. l3.5l is not satisfied. In particular, 
at the RK point, for momentum q ~ (7r, 7r, tt) + k, it has 
the form of a transverse projector, as one can see by 
Fourier transforming Eq. 13.11 



c(k) cx 



k y 



jit -y 



k 2 



k 2 



(3.8) 



This implies that only transverse excitations are gen- 
erated by af, a fact that traces back to the absence of 
monopoles, V ■ B = 0. The transverse nature of such 
excitations has already been noted by Hastings. 28 

It is important to note that these are the only gap- 
less excitations for the cubic lattice. By contrast, the 
square lattice also exhibits gapless excitations near (n, 0) 
and (0,7r), where there is a divergence of, respectively, 
s, x and Sy, signalling the immediate proximity of a crys- 
talline phase. This is consistent with the identification 
of the RK point as a critical point terminating a crys- 
talline phase. These excitations, which we have chris- 
tened piOns, are discussed further in Appendix [5] The 



absence of piOns in d = 3 signals that the RK point does 
not abut a crystalline phase but instead terminates an 
RVB phase. We will use this 'soft piOn theorem' again in 
the analysis of the FCC lattice. 



B. Effective field theory and RVB phase 

The above RK point results can now be used to con- 
strain a long wavelength field theory valid in the proxim- 
ity of this point. From it, we will be able to deduce the 
existence and properties of the neighboring RVB phase. 

To this end we follow Ref. and solve the constraint 
V • B = by writing B = V x A where A lives on the 
dual cubic lattice. As this representation comes with a 
local gauge invariance, we will pick the gauge V • A = 
to have a faithful representation of the dimer states^ 
If we now think of the time evolution generated by the 
kinetic energy, this will involve flipping plaquettes. It is 
straightforward to see that such a flip represents an (es- 
sentially) spatially local change in A. This fact and the 
RK point properties allow us to write down the action, 

S = J d 3 xdt[(d t A) 2 - p 2 (V x A) 2 - p 4 (V x V x A) 2 p.9) 

We can recognize Eq. 13.91 as the restriction to V • A = 
and Aq = of the manifestly space-time gauge invariant 
action 

S = J d 3 xdt[E 2 - p 2 B 2 - p 4 (V x B) 2 ] , (3.10) 

where E = d t A—VA , is the electric field and B — V x A 
as before. 

For p 2 — 0, this action reproduces the equal time (clas- 
sical) dimer correlators at the RK point. It also yields 
a transverse photon with ui ~ p^k 2 in agreement with 
the SMA dispersion. Further, we see that only gradients 
of B enter the Hamiltonian so that all flux sectors will 
yield degenerate ground states, as is the case in the ex- 
act solution. Finally, for p 2 < the system will be driven 
to maximize the average value of |B| — which is the case 
microscopically for v/t > 1. All of this confirms that 
Eq. 13. 101 is the correct long wavelength field theory with 
P2 changing sign exactly at the RK point. Readers famil- 
iar with this set of problems will recognize that this is in 
precise correspondence with the treatment of the square 
lattice dimer modeU 1,12 with the difference that in d = 2 
one solves the constraint by writing B = Vx in terms 
of the scalar height function. 

Unlike in d — 2, where the discreteness of the height 
field becomes relevant when p 2 > 0, in the present prob- 
lem there are no relevant operators beyond those listed 
in Eq. 13.101 when p 2 > 0. Hence for v/t < 1 the action is 
now that of the standard Maxwell theory and gives rise 
to a linearly dispersing transverse photon, uj ~ p 2 k, at 
long wavelengths. 

Other excitations: The (gapped) spinons are repre- 
sented by monomers which lead to a local divergence 
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V • B = ±1 on sublattice A/B. Hence they carry a mag- 
netic charge. As the Maxwell action is invariant under 
the duality E <-> B, we could also recast the spinons as 
electrically charged particles. Either way, they interact 
via an inverse square force law and are therefore decon- 
fined. At the RK point the vanishing of pi results in the 
force vanishing altogether. 

In addition to the photon, there is a topological defect 
of interest in the model with dimers alone. This is the 
electric monopole, being the magnetic monopole that is 
present in lattice U(l) gauge theories in our dual labelling 
of the fields. To see how it is constructed, consider solving 
the Gauss's law constraint V • E = that arises in the 
passage from (|3.1(J[) by writing E = V x P. As E lives 
on the dual lattice, P lives again on the direct lattice. 
[In quantizing the theory we take P and B to be the 
conjugate variables on a link.] Following the standard 
construction of the monopole, one can pick a classical P 
such that E has the form r/r 2 at long distances from a 
specified point of the dual lattice while a Dirac string of 
plaquettes ensures that V-E = holds everywhere— The 
monopole state then can be written down in the dimer 
representation as 

where |0) is again the ground state. Requiring that 
the resonance energy cost along the string does not 
diverge with system size quantizes the monopole flux. 
The monopole generalizes the construction of the vor- 
tex at the RK point in d = 2£i and is gapped as befits 
a deconfined phase. It was also noted in that analy- 
sis that the binding of vortices to spinons could alter 
their statistics^ Likewise in d = 3 the binding of elec- 
tric monopoles to bosonic spinons will lead to fermionic 
dyons. The monopoles also interact via an inverse square 
force. 



C. Fractionalization and quantum order 

RVB liquid phases are fractionalized in that they sup- 
port 5=1/2 spinons. In cases where there are no gapless 
excitations in the spectrum, e.g. the RVB phase on the 
triangular lattice and the quantum Hall liquids, the frac- 
tionalization is of a piece with a ground state degener- 
acy sensitive to ground state topology and a topological 
interaction between various gapped excitations — a com- 
plex of properties named topological order by Wen— In 
its pristine form, topological order is unaccompanied by 
any symmetry breaking. Most succintly, the low energy 
theory for topologically ordered systems is a topological 
field theory, the Chern-Simons theory for the quantum 
Hall states^ and the BF theory for the Z 2 RVB stated 
We will see in the next section that the FCC RVB phase 
is indeed topologically ordered in this sense. 

However, the situation with regard to the U(l) phase 
discussed in this section is different. The low energy 



theory is the Maxwell theory, which is not topological. 
Concomitantly, the gapless photon makes the distinction 
between ground states and excited states fuzzy. While at 
the RK point there is an exact degeneracy between dif- 
ferent magnetic flux sectors, elsewhere in the RVB phase, 
this degeneracy is lifted by the B 2 term in the action, the 
coefficient pi of which no longer vanishes. The splitting 
thus generated is L 3 x $ b2 /L 4 = <I? b2 /L, which van- 
ishes - only algebraically - in the thermodynamic limit 
provided $ B = o(VL). For <I? B = 0(1), the splitting has 
the same magnitude as the excitation energy of a pho- 
ton; this is no coincidence if one thinks of a photon as 
imposing a flux modulated at the photon's wavelength, 
so that locally a sector with a photon present looks like 
it belonged to a sector with a slightly different <I> B . 

As a matter of principle, one can identify the mini- 
mum energy states in the nonzero 3> B sectors as degen- 
erate ground states while reserving the term excitation 
for excited states in a given sector. Indeed, at this level 
the symmetry between E and B implies that there is an- 
other set of states with a net electric flux that also have 
an energy of 0(1/ L) and should be identified as members 
of the ground state multiplet. These are the analogs, in 
d = 3 of the height shift mode discussed by Henley for 
the RK pointiii - in d = 2 the time derivative of h is 
the electric field. Nevertheless, such degeneracies would 
be hard to detect against the background of excitations 
with similar energies and so their utility appears to be 
limited. 

Consequently, we follow Wen in identifying the 17(1) 
RVB liquid as an instance of "quantum order"— wherein 
the masslessness of the photon is attributed to the rigid- 
ity of the low energy gauge structure. 



IV. FACE-CENTRED CUBIC LATTICE: AN 
RVB Z 2 LIQUID 

The behavior of the FCC lattice quantum dimer model 
parallels that of the triangular lattice problem in d = 2. 
Like the triangular lattice, this non-bipartite lattice does 
not admit a representation by solenoidal fields. Instead, 
it exhibits topological sectors, not connected under a lo- 
cal dynamics, that are labelled by an Ising variable for 
each periodic direction on the lattice. This we call a 
winding parity to distinguish it from the U(l) winding 
numbers 3> B . The winding parity is defined in an analo- 
gous manner, by counting whether the number of dimers 
cutting the planes Sj is even or odd. 

To the right of the RK point, there exists the stan- 
dard non-flippable phase (Fig.^); a non-flippable config- 
uration can again be obtained by appropriately stacking 
staggered square lattice configurations. The identifica- 
tion of the phase for large negative v/t remains an open 
problem. 
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A. RK Point and RVB Phase 

For the RK point on the 3-torus there are 8 ex- 
actly degenerate ground states which exhibit exactly the 
same equal time correlations in the thermodynamic limit. 
More generally, there are 2 P degenerate states for a sys- 
tem with p independent non-contractible loops. The clas- 
sical analysis in Ref. LUJ has shown that the correlations 
decay exponentially to zero, with an extremely short cor- 
relation length of £ = 0.4 nearest neighbour distances. 

Turning now to collective excitations within the SMA, 
we note that there are conserved quantities, the number 
of dimers in the xy planes at wavevector q = (0, 0, tt) 
and its symmetry related counterparts. Nevertheless, 
the structure factor vanishes as well at these wavevec- 
tors, quadratically by the analyticity in momentum space 
guaranteed by the short range of the correlations in real 
space. 

To see this, consider grouping all the sites of a finite 
FCC lattice with periodic boundary conditions in the z- 
direction according to their z-coordinate. Let N z be the 
number of sites with a given z-coordinate. Now, for a 
given dimer configuration, define the number of dimers 
linking a pair of sites with neighbouring z-coordinates by 
A z z+ i. Then the number of dimers in the xy plane, n, 
at coordinate z is simply 2n(z) — N z — A ZjZ +i — A z _i ;Z . 
Therefore, for any k z = q z — tt ^ 0, 

n{q z ) ~ - ^ A Z;Z+ i exp (iq z z) [2 exp (iq z /2) cos q z ] 

z 

~ i(q z - n)A(q z ) , (4.1) 

where the last term in parentheses was expanded in k z ; 
at q z = 7r, it vanishes linearly, yielding a quadratically 
vanishing structure factor. 

From the absence of any significant crystalline corre- 
lations or gapless excitations at the RK point (no soft 
piOns) it is safe to conclude that for v/t < 1 there is a fi- 
nite RVB phase with liquid dimer correlations, an 8-fold 
ground state degeneracy on a 3-torus and a gap to all 
excitations. Indeed the first and third pieces are mutu- 
ally consistent. The vanishing of the oscillator strength 
at qo holds everywhere, while the structure factor is also 
zero exactly at qo- If the correlations remain liquid, the 
structure factor will vanish analytically at qo and hence 
the gap will persist. Evidently, matters will be different 
when a solid phase is reached. 



discussed above and the tunneling of the spinons and vor- 
tex lines lifts the ground state degeneracy by an amount 
of 0(e~ L ) and <3(e~ L ) for a system of linear dimension 
L. 



V. SUMMARY 

In this paper, we have established the existence of two 
types of RVB liquids in d = 3, which occur on the simple 
and face-centred cubic lattices. As in d = 2, the possible 
RVB phases in quantum dimer models correspond to the 
known deconfined phases in compact gauge theories and 
exhibit the same topological scaling limits. In addition 
to their direct physical interest, they allow attractively 
simple illustrations of concepts such as topological and 
quantum order. 

Remarkably, the program outlined by Rokhsar and 
Kivelson in their quest for the square lattice RVB liquid 4 
can be carried through for the simple cubic lattice in toto. 
The ability to perform calculations in a classical frame- 
work at the Rokhsar-Kivelson point to access properties 
of both ground and excited states has been invaluable for 
our study. 

Finally, it will be interesting to see if such phases can 
be realised, for example by destabilising Heisenberg Neel 
states on both lattices using frustrating exchange inter- 
actions. 



Acknowledgements 

We are grateful to Kirill Shtengel for asking stimulat- 
ing questions about the SMA analysis of quantum dimer 
models and specifically for sharing his own observations 
on the correlations of the dimer model on the square lat- 
tice. We also thank David Huse and Werner Krauth for 
collaboration on closely related work. RM would like to 
thank Leon Balents, Matthew Fisher and Mike Hermele 
for several useful discussions and SLS likewise thanks 
Vadim Oganesyan. RM was in part supported by the 
Ministere de la Recherche et des Nouvelles Technologies 
with an ACI grant, and thanks the KITP and the MPI- 
PKS Dresden for hospitality during part of this work. 
This work was supported by the NSF with grants PHY99- 
07949 (KITP), DMR-9978074 and 0213706 (SLS), and by 
the David and Lucile Packard Foundation (SLS). 



B. Topological order 

The excitations of the FCC RVB phase include the 
ubiquitous spinons and gapped vortex (vison) loops. The 
topological interaction between them is described by the 
3+1 dimensional BF action which encodes the phase fac- 
tor of e l7rni for a spinon trajectory that links n; times 
with the vortex loop>i£ Quantization of the BF action 
on closed manifolds recovers the ground state degeneracy 



APPENDIX A: CALCULATION OF OSCILLATOR 
STRENGTH 

To calculate the oscillator strength, /(q) = 
[<J%{— q), [Hqdm, crf(q)]], we need to consider the com- 
mutation properties of the dimer density operator, 
erS(q) = Y^ r crf(r) exp(iq • r) with the quantum dimer 
Hamiltonian. To make contact with the original quan- 
tum dimer paper, where the case of excitations near 
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Q = (71", tt) on the square lattice were considered^ we 
write q = Q + k; more generally, 2Q needs to be a recip- 
rocal lattice vector. 

The dimcr density operator commutes with the poten- 
tial term of Hqjjm, so that we only consider its commu- 
tators with the kinetic term. All that needs to be done 
here is to use 

[a ± ,a x ] = t2ct ± (Al) 

repeatedly. This form implies that the result in an 
isotropic phase will always be the kinetic energy operator 
(and, at the RK point, the expectation value of the loop 
fllppability) times a function of q, the details of which 
we are interested in. 

The contribution of a given square plaquette with 
dimers in the f directions at locations R and R' = R + r 
is given by 

[5£(-q) , 5?(q)]] = 

8t f exp(2iQ • R)(l + cos((Q + k) ■ r)(A.2) 

Summing this expression over all plaquettes forces 2Q to 
be a reciprocal lattice vector, and yields the expression 
proportional to the kinetic energy operator. 

The dependence on k and Q now follows straightfor- 
wardly. For example, for the case of resonons on the 
square lattice, Q = (tt, tt), and r = ±y for f = x. 
Thence, /(k) cx 1 — cosk y ~ fc 2 ~ (kxf) 2 , the original re- 
sult of Rokhsar and Kivelsonil This form also obtains for 
Q = (0, 7r), whereas near Q = (tt, 0), /(k) cx cos 2 (k y /2) 
is nonzero everywhere. 

The same analysis carries over directly to other lat- 
tices. On the cubic lattice, there are two transverse di- 
rections, and /(k) ~ (k x f ) 2 continues to hold near 

Q = (?T, 7T,7t). 

For the hexagonal case, where an elementary resonance 
move involves a triplet of dimers, one finds the resonon to 
occur near zero wavevector, Q = 0, with the same form 
of /(k)~ (kxf) 2 . 

On the triangular lattice, the same algebra yields a 
qualitatively different result. If we consider the trian- 
gular lattice as a square lattice decorated with diago- 
nal bonds pointing from the botton left to the top right 
corner of each square, the transverse projector form for 
Q = (0, 7r) is replaced by /(k) ~ fc 2 + (k y + k x ) 2 , which 
vanishes quadratically as k — ► but is nonzero in all 
directions away from k = 0. 

Finally, for the FCC lattice, where the conservation 
law involves two different dimer directions, the algebra is 
more complicated but one again finds a vanishing f(k) 
as the wavevector approaches Q = (0, 0, tt). 

APPENDIX B: SMA RESULTS IN d = 2 
1. Resonons and piOns on the square lattice 

For the bipartite square lattice, the dimer density 
cr? (q) is a conserved quantity at qo(f , q?) = z x f + qf f 




(7T,7T) 

FIG. 2: Structure factor, s^(q), for the square lattice. The 
values of q are indicated on the plot. For q y ir, the structure 
factor vanishes on the interval q x €E ] — tt, tt[. It is positive 
elsewhere, and at q = (tt, 0), it diverges logarithmically. 

for any value of qf. Thus, /(qo) vanishes. As in the cubic 
lattice, this does not lead to an entire line of zero-energy 
excitations along the Brioullin zone edge. Rather, the 
resonon only exists as a transverse excitation near (tt, tt). 

To demonstrate this, we have explicitly determined 
the structure factor for cr?(q) an d r = x. The calcu- 
lation of the structure factor for two-dimensional models 
at the RK point is straightforwardly achieved using the 
Fermionic path integrals developed for the study classical 
dimer models. A detailed calculation of correla tions f or 
the triangular lattice has been presented in Ref . 1.3-313 ll 

The result is plotted in Fig. |21 The structure factor 
vanishes everywhere on the line (q$ , tt) except for at Q = 
(tt, tt). This implies that (q, f |q, f ) = on this line away 
from (tt, tt), so that it is not a candidate for a gapless 
excitation. 

In fact, the form of the structure factor near Q is well- 
described by the transverse projector 



where q = Q + k. Together with^ /(k) ~ 1 — cosfcj, ~ 
k 2 = (kxf) 2 , this implies that there is only a single 
gapless excitation per dimer direction, namely at Q — 
(tt, tt), which is of a transverse nature. 

In Fig. [2a peak in s^(q) is visible at q! = (tt,Q). Its 
scaling with system size can be determined easily since 
the long-distance part of the correlations of the classical 
dimer model in real space are known exactly: 11-35 

CM ~ {-iy+yt^- + (-rf I . (B2) 

This form implies that the peak height grows logarithmi- 
cally with system size, and equivalently that its height a 
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distance k = q — qi from the centre decays (isotropically) 
as — log k. 

Near qi = (w, 0), the oscillator strength depends on 
k simply as /(k) ~ cos 2 (/c y /2), which is a constant to 
leading order. Therefore, E(q\,x) — 0, and there ex- 
ist further gapless excitations for each dimer direction, 
which - due to their location in reciprocal space - we 
call piOns. The dispersion obtained within the SMA is 
logarithmic, but more on that below. 

We mention in passing that the same analysis carries 
over to the hexagonal lattice, where the resonon exists at 
Q = (0, 0) and the now inappropriately named piOn at 
(47r/3,0). 

This similarity is efficiently encoded in the language 
of two-dimensional height models, analysed in detail by 
Henley^i Both resonon and piOn can be connected back 
to slow modes in the height model at zero wavevector, the 
connection to the dimer modes being via the gradient and 
vertex operators, respectively. Indeed from the height 
analysis one learns^ that the structure factor at qi + k 



decays in imaginary time as e~ fcv ^ which confirms that 
there are gapless excitations near qi but also shows that 
the SMA does not do a very good job of constructing 
them. 



2. Absence of gapless modes on the triangular 
lattice 

The argument that the SMA mode for on the trian- 
gular lattice is gapped at all momenta, despite the pres- 
ence of a conservation law, is completely analogous to 
the one presented above for the FCC lattice. It again 
involves showing that the corresponding structure factor 
vanishes at least quadratically. The relevant momentum 
for dimers pointing in the x direction is Q = (0, tt). A 
direct computation of the structure factor using the re- 
sults of Ref . l33l yields the same result. This SMA result 
was already noted in Ref. 
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